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1.  INTRODUCTION 


Significant  advancements  in  power  electronic  devices,  electric  machine  design,  and 
power  converter  packaging  have  made  it  possible  for  the  U.S.  Navy  to  consider  electric 
drive  for  its  new  land-attack  destroyer,  DD-21.  Many  of  the  benefits  of  electric  drive  stem 
from  incorporating  it  into  an  Integrated  Power  System  (IPS),  where  both  ship  service  and 
ship  propulsion  power  are  derived  from  a  common  set  of  prime  movers.  The  commercial 
marine  industry,  including  Alstom,  ABB  and  Siemens,  has  exploited  the  benefits  of  IPS  for 
more  than  a  decade  in  cruise  liners,  ferries,  tankers,  and  research  vessels.  The  advantages 
of  IPS  are  manifold  including;  reduced  lifecycle  costs,  flexible  ship  layout,  cross-connect 
capability,  reduced  signature,  optimized  prime  mover  operation,  simple  and  rapid  propeller 
control,  and  a  pool  of  available  electric  power  for  future  electric  auxiliaries  and  weapons. 

In  addition  to  the  innovation  of  IPS,  the  U.S.  Navy  is  also  considering  new  paradigms 
in  power  distribution.  In  particular,  traditional  AC  radial  distribution  architectures  are 
being  replaced  by  AC  zonal  architectures.  The  concept  of  a  DC  Zonal  Electric 
Distribution  System  (DC  ZEDS)  is  being  investigated  and  advocated  for  its  improved 
fault-response  speed  and  distributed  intelligence.  Initiatives  such  as  electric  drive,  IPS  and 
DC  ZEDS  depend  on  the  maturation  of  certain  enabling  technologies.  A  partial  list 
includes  power  devices,  compact  high-bandwidth  power  converter  topologies,  high-power 
low-speed  electric  machinery,  reliable  automation  technology,  and  advanced  control  and 
decision-making  algorithms. 

The  research  effort  reported  on  in  this  document  concentrates  on  developing  models 
and  analyses  required  to  control  multiple  power  inverters  operating  in  parallel.  One 
potential  application  for  such  work  is  the  main  Propulsion  Motor  Module  (PMM).  The 
PMM  may  consist  of  an  AC  motor  driven  by  a  DC- AC  power  inverter.  One  option  for 
such  an  inverter  is  to  employ  a  high-power  slower-switching  main  inverter  in  parallel  with 
a  lower-power  fast-switching  inverter.  The  rationale  for  such  a  choice  is  explained  as 
follows.  A  typical  surface  combatant  may  require  on  the  order  of  100,000  SHP  to  achieve 
full  speed  (50,000  SHP  per  shaft),  indicating  a  PMM  power  rating  of  approximately 
40MW.  In  order  to  achieve  high-fidelity  machine  waveforms  and  rapid  control,  high- 
bandwidth  (fast-switching)  power  converters  are  required.  At  high  power,  this  implies 
IGBT  technology.  Unfortunately,  to  achieve  20kHz  or  greater  device  switching,  IGBT- 
based  inverters  are  limited  to  about  2MW.  Despite  the  fact  that  the  machine  may  be 
divided  into  more  than  three  separate  windings  (phases),  this  would  still  require  multiple 
paralleled  modules  and  an  unacceptably  large-volume  solution.  For  high  power,  devices 
must  then  be  series  connected,  resulting  in  much  lower  admissible  switching  frequencies. 
One  compromise  then  is  to  provide  the  "bulk"  power  of  the  drive  using  a  slow-switching 
high-power  inverter.  The  waveform  fidelity  could  then  be  assured  by  paralleling  with  a 
fast-switching  lower-power  unit  which  may  have  to  process  no  more  than  10%  of  the  rated 
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power.  Thus  the  main  "bulk"  inverter  would  be  rated  at  ~36MW  and  the  "high-fidelity" 
inverter  would  be  rated  at  “4MW,  theoretically  realizable  by  two  units  in  parallel.  The 
"high-fidelity"  unit  would  "fill  in"  the  waveform  distortion  produced  by  the  "bulk"  unit 
which  is  constrained  to  switch  at  l-2kHz. 

Another  application  that  may  require  paralleled  inverters  is  DC  ZEDS.  In  DC  ZEDS, 
DC-to-AC  inverters  convert  the  interzonal  DC  voltage  to  the  necessary  AC  voltage 
amplitudes  and  frequencies.  To  avoid  developing  and  inventorying  inverters  of  many 
different  power  ratings  and  to  implement  redundancy  in  the  distribution  system,  inverters 
are  paralleled.  High-switching  frequency  units  are  required  to  achieve  suitable  AC 
waveforms  to  optimize  efficiency  and  minimize  the  acoustic  signature  of  auxiliary  systems. 

In  order  to  tackle  the  problem  paralleling,  several  questions  must  be  posed:  Will  the 
units  be  co-located  and,  therefore,  can  they  share  information?  If  they  are  completely 
autonomous,  how  are  they  synchronized  and  how  are  the  reference  settings  (i.e.,  voltage, 
current)  determined?  What  strategy  will  be  used  to  ensure  proper  load  sharing?  The 
following  effort  imposes  the  following  boundaries  on  this  problem:  First,  the  units  are 
assumed  to  be  co-located  as  illustrated  in  Figure  1.  Note  that  in  Figure  1  the  three-phase 
outputs  of  each  unit  are  designated  by  a  single  connection  line.  Coupling  inductors  are 
used  to  aid  in  proper  sharing  and  to  limit  unwanted  circulating  currents  between  the  units. 
A  common  DC  voltage  input  is  then  applied  to  each  individual  inverter  unit.  Second,  an 
inner  current  control  is  assumed  regardless  of  the  type  of  load.  With  identical  inverter 
units,  this  implies  that  each  unit  should  produce  an  equal  portion  of  the  commanded  load 
current.  Third,  a  master  control,  also  co-located,  will  produce  the  commanded  current 
signals  for  each  unit  and  establish  the  common  electrical  angle  to  ensure  synchronization. 
As  a  consequence,  if  a  fault  occurs  at  one  unit  or  it  fails,  the  master  control  is  responsible 
for  re-allocating  the  current  requirement  between  the  remaining  units  within  the  bounds  of 
their  ratings.  Furthermore,  by  assuming  co-located  units,  the  assertion  of  a  co-located 
master  controller  is  not  unreasonable.  Also,  this  does  not  prohibit  a  fail-back  local  droop 
control  should  the  master  fail.  The  assumed  "global"  or  master  control  is  indicated  in 
Figure  2,  where  separate  q-  and  d-commanded  currents  are  specified  in  the  synchronous 
reference  frame  (outlined  later  in  this  document).  The  quantity  co*  is  the  commanded 
electrical  radian  frequency  which  could  be  a  constant  or  could  be  a  variable,  say  in  a 
variable-speed  drive.  The  fourth  and  final  assumption  is  that  the  local  control  will 
appropriately  handle  circulating  or  zero-sequence  currents.  The  local  control  refers  to  how 
the  commanded  currents  at  each  "local"  unit  are  converted  into  actual  switch  gate  signals. 
One  approach  to  such  a  local  control  is  illustrated  in  Figure  3.  The  actual  inverter  output 
currents  (iai,  ibi,  id)  are  fed  back  to  the  local  control.  After  being  transformed  into  the 
synchronous  reference  frame,  they  are  compared  with  the  desired  current  values  ("desired" 
is  designated  by  an  asterisk). 
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Inverter  Module 


Figure  1.  Inverter  module  containing  five  inverter  submodules 


Figure  2.  Synchronizing  and  sharing  global  control 
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Figure  3.  QDO  local  control 


The  qdO-controllers  produce  modulating  signals  which  are  then  inverse  transformed  back 
into  "physical"  (abc)  variables  which,  for  instance,  can  then  be  used  in  a  standard  Sine- 
Triangle  Pulse- Width-Modulation  (STPWM)  algorithm.  Other  approaches  including  Space 
Vector  Modulation  are  easily  substituted  into  the  block  diagram. 

The  principal  goal  of  this  work  is  to  establish  a  theoretical  basis  for  the  stability  of  a 
particular  control  implementation,  while  in  the  process  evolving  the  necessary  system 
modeling  equations  required  for  detailed  simulation.  The  initial  development  is  conducted 
for  a  generic  three-phase  Resistor-Inductor  (R-L)  load,  but  since  many  different  three-phase 
loads  are  admissible  for  shipboard  applications,  the  results  are  then  extended  to  a  generic 
three-phase  induction  machine  and  a  generic  three-phase  permanent-magnet  synchronous 
machine.  The  following  details  are  presented  in  this  document: 

•  the  dynamic  equations  are  developed  for  two  inverters  powering  an  R-L  load 
(given  an  assumed  local  current  control); 

•  an  equivalent  single-inverter  system  is  derived  and  gains  are  analytically 
determined; 

•  a  mechanism  for  mapping  the  single-inverter  gains  to  a  multiple  number  of 
paralleled  inverters  is  outlined  and  verified; 

•  the  local  zero-sequence  control  is  designed  to  handle  components  of  the 
potential  circulating  currents; 

•  the  above  results  are  extended  to  a  symmetrical  three-phase  induction  machine 
load; 


-  5  - 

•  the  above  results  are  extended  to  a  round-rotor  three-phase  permanent-magnet 
synchronous  machine  load; 

•  the  necessary  equations  for  the  vector-control  of  induction  machines  and 
permanent-magnet  synchronous  machines  are  outlined; 

•  the  design  of  an  outer  speed-control  loop  appropriate  to  each  machine  type  is 
presented. 

Much  of  the  aforementioned  development  is  equation  and  reference  frame  intensive.  Many 
intermediate  steps  are  included  to  provide  guidance  for  any  future  modifications  or 
extensions  of  this  work  (i.e.,  different  loads  or  different  control  strategies). 
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11.  TWO  PARALLELED  INVERTERS  POWERING  AN  R-L  LOAD 

The  basic  system  under  consideration  is  illustrated  in  Figure  4.  Note,  the  inverters 
share  a  common  DC  source  bus  and  each  is  coupled  to  the  R-L  load  through  inductors 
labeled  Li  and  L2.  The  voltage  ,  is  referenced  from  the  center  of  the  a-phase  leg  (of  unit 
1)  to  a  fictitious  mid-point  of  the  DC  voltage  supply  (designated  o).  These  voltages  are 
sometimes  referred  to  as  "pole  voltages."  The  load  is  assumed  to  be  wye-connected  with 
the  neutral  point  (labelled  n)  left  floating. 


Figure  4.  Two  paralleled  inverters  powering  an  R-L  load 
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Since  the  interconnection  of  inductances  L|,  Lz  and  Ll  form  a  cutset,  all  three  inductor 
currents  cannot  be  state  variables.  Thus,  one  is  eliminated  by  virtue  of  the  algebraic 
constraint.  If  we  choose  the  currents  in  L]  and  Lz  as  our  state  variables,  the  dynamics  of 
the  circuit  portrayed  in  Figure  4  can  be  established. 


A.  Circuit  Representation  in  Physical  Variables 

The  following  Kirchhoff  s  Voltage  Law  (KVL)  loop  equations  may  be  written  for  the 
network: 


“Vao,l  +  +  RL(ial+ia2)  +  Ll” (ial  +  ia2)  +  Vno  -  0 

d  d 

“Vbo,l  +  ^L(ibl+ib2)  +  Ll  ( ibl  +  ib2 )  +  Vno  =  0 

d  d 

“''col  +  Lj— id  +  RL(ici  +  ic2)  +  Ll  — (ici+ic2)  +  Vno  =  0 
d  d 

“Vao,l  +  ^  ~  ^ 

T  d  •  T  d  . 

“■Vbo,!  +  “  L2  — lb2  +  Vbo,2  -  0 

“Vco.I  +  Ll-— id  “  ^^L2  'Vco,2  =  0 

Equations  (l)-(6)  may  be  expressed  in  vector  notation  as: 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


V  d  “7^  V  w  ,  d  “1:4 

“Vabco,I  +  M— labc,I  +  ^LCJabc,!  +  labc,2)  +  J-L(‘^labc,l 


dt 


labc,2)  + 


=  d 


(7) 


-y. 


abco,l 


+  Ll 


d  ^ 
dt 


L2 


d 

dt 


^abco,2  “  ^ 


(8) 


where  for  instance  Tabc,i  =  [iai.ibi.id  f  and  Vabco,i  =  [Vao.i.Vbo,i.Vco,i  The  voltage  v„o  is  the 
voltage  from  the  neutral  point  of  the  load  back  around  to  the  mid-point  of  the  DC  voltage 
supply. 


B.  Circuit  Representation  in  the  Synchronous  Reference  Frame 

These  equations  may  be  transformed  into  the  synchronous  reference  frame  (designated 
by  the  superscript  "e")  using  the  identities 


^qdO,l  “  Rs  ^abco.l 
'^qd0,2  ~  Rs  ^abco,2 


iqdO.l  “  iflbc,! 
iqd0,2  “  Rs  iabc,2 


(9) 

(10) 

(11) 

(12) 
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where  the  diffeomorphic  transformation  matrix  Kf  is  given  by 


= 


2  »  2  2n  .  2  /o  ,271  . 

jCOsGe  jCOSCGe-  — )  J  005(05+  — ) 

2-0  2  .  27t  2  -  2jc 

ySinGe  jsm(0e-  — )  jsin(0e+  — ) 


(13) 


and,  for  instance,  the  notation  implies  that  ^qdoj  =  [Vq®|,Vd®,,vo.i  ]^.  The  zero-sequence 
quantity  (i.e.,  vq,,)  is  generally  not  given  a  superscript  since  it  is  equivalent  in  every 
reference  frame.  The  angle  0e  varies  at  a  radian  frequency  of  cOg,  which  for  instance,  would 
be  the  fundamental  or  desired  output  frequency  of  the  inverters.  Solving  for  the  "abc" 
quantities  in  (9)-(12)  and  then  substituting  into  (7)-(8)  gives: 


-(Ks'^r'^qVl  +  (L,+Ll)  — 


dt 


(IQ'r'iqdO,.  +  RL(Ks')-‘iq=d0,,  + 


ev-iye 


Rl  ( )■  %0.2  +  Ll  -^  [  ( )-  %o.: 


-(K/)-'^qVi  +  (K3"r'iqV 


(K3*)-'?q'd0,2  +  ^ 


Multiplying  (14)“(15)  by  and  employing  the  product  rule  illustrated  in  (16) 


dt 


C,V.]  =  (Ks'^r'-^tfdo.i  + 


(14) 


(15) 


(16) 


gives 


-VqWl  +  (L,+LL)-|-TqVl  +  (L,+LL)K/-|-[(K3^r']tqVl  + 


Rl  (  iqdO,  1  +  iqdO.2  )  +Ll  Kj®  —  |^  (  K,®  )  ' 


*qd0,2 


=  '6 


-^qdO.I  +  L,-|-tqV,  +  L,K;-|- 


ex-1 


(K/y 


lqdO,I  ” 

'■^[(Ks*^)  'j  jqV2  +  '^qdO.2 

Equations  (17)  and  (18)  are  simplified  by  observing  that 


d  -T>e 

-2  lqdO.2 


re  ,  ^ 


Ks^-^(K3^r'  = 


0  (Oe  0 

-cOg  0  0 
0  0  0 


(17) 


(18) 


(19) 
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and 


K® 


- 1 

< 

■  o' 

Vno 

il 

0 

< 

.V„o 

(20) 


Thus,  expanding  the  resultant  vector  equations  gives  the  representation  of  the  two-inverter 
system  in  the  synchronous  reference  frame: 


“Vq®]  +  (Li+Ll)  — iq®l  +  (Lj  +  LL)0)eidj  + 


( iqtl  +  iq!2 )  +  Ll  —  iq®2  +  —  0 


■V<fl  +  (Ll+Lbl'^idJ  ~  (Ll  +  LL)®eiq!l  + 


RL(id!i  +  ida)  +  Lf^idta  “  Lt^eiqta  -  0 


•^0,1  +  (Li+Ll)  — io,i  Rl ( io,i  +  io,2 )  +  LL'^io,2  +  =  0 


^  •e  .T-.  r 


-Vq!l  +  Ll— iqj  +  L,C0eid!l  -  ^2 iq!2  “  L2®eid!2  +  Vq®2  =  0 
"''d!:  +  “  l-l^eiq!]  ~  L2  id!2  +  +  ^d,2  =  9 

T  d  .  ,  d  . 

"Vo,!  +  L]  —  lo,i  -  L2  — lo,2  +  Vo, 2  -  0 


(21) 


(22) 


(23) 

(24) 

(25) 

(26) 


Recall,  the  voltages  Vq®i  and  vj®,  are  the  transformed  inverter  #1  pole  voltages  while  vo,i  is 
the  zero-sequence  voltage  for  that  unit. 

With  a  wye-connected  load,  Kirchhoff  s  Current  Law  (KCL)  demands  that 


( ial  +  ia2 )  +  ( ibi  +  ib2 )  +  ( id  +  ic2  )  =  0  (27) 

But  from  the  identities  (11)-(12),  the  inverter  zero-sequence  currents  relate  to  the  output 


currents  by 

i  1 

«o..  -  j  1 

ial  +  ibI  +  id  j 

(28) 

1  1 

'' 

(29) 

10,2  -  J  1 

ia2  +  ib2  +  ic2 

SO  that  (27)  then  implies  that 

io,2  =  ~io,i 

(30) 

Therefore,  the  inverter  #2  zero-sequence  current  is  NOT  an  independent  state  variable,  but 
instead  is  algebraically  related  to  io  i.  This  is  a  direct  consequence  of  the  wye-connection 
of  the  load  with  the  neutral  floating.  Substituting  (30)  into  (26)  gives 
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-Vo.i 


+ 


+ 


+  ''0,2 


0 


which  simplifies  down  to 


(31) 


d  .  _  VqJ  ~  ^0,2 

dt  “  L,  +  La 


(32) 


Thus,  any  difference  between  the  developed  zero-sequence  voltages  of  the  units  will  result 
in  a  zero-sequence  current,  which  must  flow  between  the  units.  Furthermore  substituting 
(30)  into  (23)  results  in 

~''0,l  +  +  ''no  =  0  (33) 


which  allows  us  to  solve  for  v^o  by  first  substituting  (32) 


Vno  =  Vo,l 


L, 

Li  +  La 


V0,I  -  Vo,2j 


(34) 


or  more  compactly 


V 


no 


L2''0,1  +  Ll ''0.2 

Li  +  La 


(35) 


The  dynamic  governing  the  zero-sequence  current  is  described  by  (32)  while  the 
remaining  state  variable  equations,  (21)-(22)  and  (24)-(25),  may  be  rewritten  in 
vector/matrix  form  as 


and 


La~iq!l2  -  RaiqJ2  ~  ®eLaidJ2  +  Q'^qJ2 


dt 

dt 


La~id*12  -  ®eLaiqJ2  +  Raid*I2  +  Q^dJ2 


where 


'q,12  - 


’q*i 

Liq'aJ 


Id, 12  - 


id'i 

Lid!2j 


^q^l2  = 


Vq!2j 


Vd,12  = 


u  = 


L|  +  Ll  Ll 

L,  -La 


(36) 

(37) 

(38) 

(39) 

(40) 

(41) 

(42) 
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Ra  = 


(43) 


Ca  = 


1  0 
1  -1 


(44) 


Equations  (36)-(37)  may  be  placed  into  state-space  normal  form 

■^iq!l2  =  La'Raiq!l2  “  ®eid!l2  +  ^  '  Ca7q®i2 


d  ^ 


■^id'l2  =  C0eiq42  +  La’Raid!l2  +  ^ 

where  for  the  reader’s  information 
RlL2  RlL2 


^In 


La''  Ra  = 


Lr‘  c,  = 


Lden  ^den 

RlLi  RlLi 


Lden  Ljei, 

-L2  -  Ll 


Lden 

L-l 

Lden 


Lden 

-L,  -  Ll 


■"den 


and 


(45) 

(46) 


(47) 


(48) 


Laen  =  -L,L2  -  LlL2  -  L,  Ll  (49) 

Thus,  from  (45)-(46)  we  observe  that  the  "q"  and  "d"  dynamics  are  coupled  by  the  "speed 
current"  terms  and  cOeTd!^- 


C.  Synchronous  Reference  Frame  Control  Equations 

In  order  to  investigate  the  stability  of  a  proposed  control  architecture  and  to  select 
suitable  control  gains,  the  control  dynamics  must  be  modeled  and  linked  to  equations  (45)- 
(46).  Assume  that  each  inverter  unit  is  supplied  a  common  synchronous  reference  frame 
angle  0e,  which  can  then  be  used  locally  together  with  the  inverter  output  currents  to 
calculate  the  inverter’s  qd-currents  (see  Figures  2  and  3).  These  currents  may  then  be 
compared  against  desired  qd-synchronous  reference  frame  currents  and  processed  through  a 
Proportional-plus-integral  (PI)  controller.  The  outputs,  still  mathematically  in  the 
synchronous  reference  frame,  would  then  be  inverse  transformed  back  into  abc-frame 
(physical)  variables,  which  could  serve  as  the  modulating  signals  in  a  Sine-Triangle  Pulse- 
Width-Modulation  (STPWM)  strategy  (see  Figure  3).  A  Pl-control  is  used  since  for 
sinusoidal  reference  abc-currents,  the  commanded  qd-currents  are  constants  in  the 
synchronous  reference  frame.  Thus,  the  integral  action  will  guarantee  zero  steady-state 
error,  while  the  proportional  action  guarantees  improved  stability  margin.  The 
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implementation  of  the  integrator  must  account  for  "integrator  windup"  which  may  occur 
when  large  operating  point  changes  are  specified. 

Given  a  bipolar  triangle  waveform  varying  between  A,  and  -A„  the  control  is  said  to 
remain  in  the  "linear  modulation"  range  as  long  as  each  modulation  signal  has  magnitude 
less  than  A,.  If  the  control  is  constrained  to  operate  in  the  "linear  modulation"  range,  then 
the  following  relates  the  fundamental  modulation  signal  to  the  respective  fundamental 
inverter  output  line-to-neutral  voltage  (as  long  as  the  fundamental  frequency  is  much  lower 
than  the  carrier  frequency): 


Val,n 

Vbl,n 

^amod.l 

"^dc 

“  2A, 

'^bmod.  1 

Vcl.n 

^cmod.  1 

(50) 


The  frequency  spectra  for  the  line-to-neutral  voltages  is  identical  to  the  pole  voltages 
except  for  the  triplen  harmonics.  Therefore,  the  fundamentals  are  identical  and  would  thus 
be  also  related  by  (50).  Transforming  (50)  into  the  synchronous  reference  frame  gives  us 
the  "average"  qd-voltage  relationships: 


e 

Vdc 

'^qmod,  1 

e 

Vdl 

2  At 

^dmod.I 

Voj_ 

^Omod.l 

(51) 


Vdc 

The  term  Kp^m  =  is  termed  the  STPWM  gain  and  is  clearly  a  function  of  the  inverter 

input  DC  voltage  and  the  amplitude  of  the  carrier  triangle  waveform  (clearly  this 
modulation  may  take  place  via  a  Digital  Signal  Processing  implementation  and  therefore 
the  amplitude  A,  is  dictated  by  the  programmer).  The  dynamics  of  the  Pl-control  and  the 
interface  with  the  modeling  equations  (45)-(46)  are  therefore  given  by: 


(52) 


;e*  :  e 

^d,l  “  Id.l 


(53) 


d 

dt 

-  iq.2  - 

dt 

=  id,2  - 

id"2 

''q^l 

= 

^pwm  ^pq 

;e*  ^ 

' 

VdJ 

= 

^pwm  ^^d  1 

^d,l 

■  idj 

= 

^pwm  ^pq  j 

’iq^j- 

“  iq^2 

^iq  ^q,  I 
^d.I 

kpwm  ^iq  ^q,2 


(54) 

(55) 

(56) 

(57) 


(58) 
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'^d,2  ~  l^pd  *d,2  id, 2  l^id  ^d, 


where  the  quantities  with  "*"  denote  commanded  or  desired  values.  The  resultant  set  of 
equations,  (45)-(46)  and  (52)-(59)  provide  the  8th-order  system  of  equations  in  state-space 
normal  form.  As  long  as  the  average  qd-voltages  adequately  represent  the  bulk  of  the 
dynamics,  the  eigenvalues  of  the  resultant  system  matrix  characterize  the  stability  of  the 
interconnected  inverters  and  load.  Reiterating,  these  equations  would  not  be  used  for 
detailed  simulation  since  we  have  average-value  modeled  through  the  inverter;  however, 
the  other  purpose  of  modeling  is  to  facilitate  control  design,  which  is  what  we  are  setting 
the  stage  to  do  here.  Those  equations  take  the  form: 


■^iq!l2  -  A-iiq®i2  -  0)eidtl2  +  KiqA2'?q®l2  +  KpqA2iqj2 


■^>d!l2  -  ®ejq!l2  +  ^11(^12  +  Kjj A2'j?d!l2  +  KpdA2id!l2 


■  e  _  -j»e  fe* 

-  ~lq,l2  Jq,I2 


Q  ^  p  *T>P  * 

=  “Id, 12  +  ld,12 


where  the  "primed"  gains  are 


^pq  ”  ^pwm  ^pq 


Kiq  —  Kjq 


^pd  ““  ^pwm  ^pd 


Kjd  —  ^id 


Technically,  as  will  be  discussed  further  in  this  document,  the  primed  gains  may  also  be 
used  to  represent  additional  scaling  effects  such  as  the  non-unity  gain  of  current  sensors. 
The  matrices  in  (60)-(61)  are  given  by 

~  R^L2.+  K;^(U  +  L0  R^L^-K^qL^ 

_  hjejj  Ljen 

=  RLL,-KpqLL  RLki+K^qCLi+Lt) 

hden 


L2  +  Ll 

Ll 

hden 

Lden 

Ll 

Li  +  Ll 

Lden 

Lden 

RlL2  +  Kpd(L2  +  Ll)  “  KpjLl 


RlL]  ~  Kp^L^  Rl L]  +  Kpd  ( L|  +  Ll  ) 

bden  Lden 
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The  notation  in  (60)-(63)  uses  the  vectors 


The  resultant  system  matrix  takes  the  form: 


A] 

~COel2x2 

Kiq  A.2 

02x2 

“ei2x2 

^2x2 

Km  A; 

~l2x2 

^2x2 

02x2 

02x2 

^2x2 

^^2x2 

02x2 

02x2 

(71) 

(72) 

(73) 

(74) 

(75) 


where  02*2  is  a  two-by-two  null  matrix  and  12^2  is  a  two-by-two  identity  matrix.  The 
eigenvalues  of  (75)  dictate  the  stability  and  strongly  influence  the  transient  response 
characteristics  of  the  system.  We  will  illustrate  in  subsequent  sections  how  manipulating 
the  Pl-control  gains  changes  the  eigenvalue  locations. 
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III.  DETERMINING  THE  CONTROL  GAINS 

In  order  to  establish  a  set  of  gains  for  each  unit,  let’s  consider  converting  the  two 
separate  inverter  units  into  a  single  quivalent.  To  do  this,  the  per-phase  steady-state 
equivalent  circuit  of  Figure  5  is  helpful.  The  Thevenin  Equivalent  circuit  looking  back 
from  the  load  is  depicted  in  Figure  6.  If  we  assume  near-identical  units  (L|=L2  and  equal 
voltages),  then  the  Thevenin  circuit  reduces  to  simply  the  inverter  voltage  in  series  with 
half  of  the  coupling  inductance. 


Rl 


V  ^2 
L1  +  L2 


L| 

L1+L2 


Figure  6.  Per-phase  steady-state  Thevenin  circuit 
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Thus,  to  represent  two  units  by  one  equivalent,  the  connection  inductance  used  for  the 
equivalent  must  be  one-half  the  inductance  used  for  either  channel.  The  following 
procedure  is  adopted: 

•  model  the  equations  for  a  single  inverter  with  connection  inductance  given  by 
Lx  =  Y  L,  (thus,  for  2  units  in  parallel); 


•  given  desired  closed-loop  pole  locations,  substitute  the  control  model  into  this 
representation  and  solve  for  the  two  q-channel  and  the  two  d-channel  gains; 

•  modify  the  gains  to  account  for  the  fact  that  two  inverters  are  used.  This 
happens  because  each  unit  is  only  responsible  for  supplying  half  of  the  desired 
q-  and  d-axis  load  current. 


To  better  illustrate  the  last  bullet,  a  single  aggregate  inverter  would  have  a  q-axis 
modulating  voltage  established  by: 

t 

V|,„,0d  =  Ki  J  (  if  -  )  dt  +  Kpq  (  i  =  ‘  )  (76) 

0 


In  a  system  with  two  paralleled  inverters,  the  commanded  current  is  dividely  evenly 
between  the  units  so  that  the  q-axis  modulating  signal  for  unit  one  would  be 


^q.modl 


t  e*  -e  |e*  :e 

=  Kjq,  +  Kpq,(-^-^) 


(77) 


or  to  make  the  point  more  explicit 


Ki, 


'^q.modl 


qi 


jo;-- 


ia")clt 


K, 


pql 


(if -ip*) 


(78) 


Thus,  for  one  of  the  two  parallel  units  (designated  with  subscript  1)  to  produce  the  same 
voltage  as  the  equivalent  unit  (no  subscript  number),  it  must  follow  then  that 


Kip 

Kpp 


Kjql 

2 

^pqi 

2 


(79) 

(80) 


Therefore,  upon  finding  the  gains  for  a  single  unit,  the  gains  for  each  of  the  two  separate 
units  are  found  by  doubling  those  values.  This  point  will  be  revisited  during  a  subsequent 
example. 


A.  Equivalent  Inverter  System  Equations 

The  development  of  the  equations  required  to  find  the  single-unit  gains  is  presented 
next.  The  inductance  Lx  represents  the  single  unit  (or  equivalent  connection)  inductance. 
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Omitting  consideration  of  the  zero-sequence  dynamics  (since  with  a  wye-connected  load 
and  a  single  inverter,  there  would  not  be  any),  the  voltage  equations  in  the  synchronous 
reference  frame  are  given  by: 

-v|  +  (L,  +  LL)-|-iq^  +  (4  +  LL)a)ei|  +  =  0  (81) 

-v|  +  (L,  +  LL)-^i|  -  (L,  +  LL)«)ei|  +  RlII  =  0  (82) 

or  in  state-space  form 


d  -e  e 

—  Iq  -  -©eld  - 

d  R|. 


U  +  L, 


1 


dt'^'  L,  +  Ll 


Id 


Lx  +  Ll 
1 


Coupling  these  with  the  Pl-control  equations 


X  =  i  —  i 
(jj  M  ‘q 


dt 


^d 


Kp„:„[Kpq(i®*-i|)  +  KiqX®] 

=  Kp,^[Kpd(id®*-i|)  +  KidX|] 
yields,  in  state-space  normal  form, 

-(Rl+k;,) 


dt 


Lx  Ll 

®e 

-1 

0 


K, 


■pq 


Lx  +  Ll 
0 


K, 


■pd 


1 

0 


Lx  +  Ll 
0 
1 


Ly  +  Ll 


•Vd 


Ki 


•iq 


-(Rb  +  Kpd) 

Lx  Ll 
0 

-1 


Lx  +  Ll 
0 


0 

0 


0 

Kid 


Lx  +  Ll 
0 
0 


where  K^q  =  Kp^^Kpq  and  Kj,  =  Kp„„Kiq  have  been  used  for  notational  simplicity. 

B.  A  qd-Current  Loop  Design  Strategy 
Equation  (89)  can  be  further  simplified  if  we  assign 


Kpq  = 


Rl  +  K, 


■pq 


(83) 

(84) 

(85) 

(86) 

(87) 

(88) 


(89) 


Lx  +  Ll 


(90) 
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K 


pd 


+  l^pd 
Lx  Ll 


,  _  Kj, 

L,  +  Ll 


Km  = 


K 


‘id 


Lx  +  Ll 

The  characteristic  equation  for  (89)  is  then  found  to  be 


x*  +  +  +  (Kp;,K;<,  +  (o2  +  Ki;  +  Kd)^^ 

+  (  Kpq  Kd  +  Kpd  Kq  )  ^  +  Kjq  Kjd  =  0 
If  this  expression  is  set  equal  to  a  desired  characteristic  equation 

+  d3A,^  +  dll'?"  +  d|  A.  +  do  =  0 


(91) 

(92) 

(93) 

(94) 

(95) 


coefficient  matching  yields  the  following  nonlinear  coupled  set  of  equations: 


Kpq  +  Kpd  =  d3  (96) 

K;qKpd  +  CD2  +  Ki;  +  Kd  =  d2  (97) 

K;qK>K;dK"  =  d,  (98) 

K"Kd  =  do  (99) 


A  Newton-Raphson  algorithm  was  selected  to  iteratively  solve  these  equations  to  find  the 
"double-prime"  unknowns  from  which  the  actual  gains  can  then  be  derived.  The  first  step 
is  to  re-express  (96)-(99)  in  terms  of  the  following  four  nonlinear  functions: 


=  Kpq  +  Kpd  -  d,  (100) 

h  =  k"  k;;  +  ©2 + K"  +  Kd  -  d2  (101) 

f,  =  K;,K>K;dKq  -  d,  (102) 

fo  =  KqKd  -  do  (103) 


The  update  equation  is  given  by 


■*^pq,new 

^pd,new 

^pq,old 

^pd,old 

5?  .X."  ^ 

-1 

fsCCid) 

f2(Cid) 

^iq.new 

^id,new 

^iq,old 

^id,old 

sr  ^ 

_ 1 

(104) 


where  1^=  [f3,f2.fi.fo]^  and  it"  =  [  K^,  K^'o,  K',  Kd  The  notation  ^’,d  implies  that  it"  is 
evaluated  at  the  old  (or  previous)  value.  The  specified  Jacobian  matrix  is  given  by 
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5^ 

5i^" 


(<ld)  = 


110  0 
l^pd.old  l^pq.old  ^  ^ 

l^id.old  ^iq,old  l^pd.old  I'^.old 
0  0  Kid. old  Kjqoid 


(105) 


C.  Design  Example 

Given  the  following  desired  Bessel  closed-loop  pole  locations: 

=  8000(-0.6573  ±  j 0.8302)  =  -5258.4  ±  j 6641.6  (106) 

A,3,4  =  8000(-0.9047  ±  j0.2711)  =  -7237.6  ±  j2168.8  (107) 

(which  provide  a  closed-loop  bandwidth  of  about  8000rad/sec  or  1273Hz),  the  desired 

coefficients  of  the  characteristic  equation  (95)  are 


dj  =  2.4992000e+04 

(108) 

d2  =  2.8108095e+08 

(109) 

d]  =  1.6391316e+12 

(110) 

do  =  4.0966232e+15 

(111) 

The  authors  selected  a  bandwidth  of  1300Hz,  so  that  the  control  would  not  interact  with  an 
assumed  switching  frequency  of  at  least  20kHz.  If  additional  signal  filtering  is  used  (i.e., 
for  the  local  current  measurements),  then  the  control  bandwidth  may  need  to  be  reduced  to 
prevent  unwanted  interaction  with  the  filters.  This  would  need  to  be  explored  via  detailed 
simulation  and  possibly  "tweaked"  during  implementation.  Update  equation  (104)  is  run  in 
MATLAB  with  (Oe= 377 rad/sec  yielding  the  following  double-primed  gains: 


Kpq  =  1.0485e+04 

(112) 

Kpd  =  1.4507e+04 

(113) 

Kiq  =  7.16866407 

(114) 

Kid  =  5.71476407 

(115) 

Note,  equation  (104)  is  not  at  all  parameter  dependent  and  is  only  a  function  of  the  form  of 
the  system  equation  (89).  The  actual  gains  for  the  single-unit  controller  are  determined 
from 


Kpq  = 


Kpd  = 


Kiq  = 


(L,  +  LL)Kpq  -  Rl 
Kpwm 

(L,  +  LL)K;d  -  Rl 


Kpwm 

(L^  +  LOKj; 
Kpwm 


(116) 

(117) 


(118) 
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Kid 


(4  +  LL)Ki'd 


K, 


pwm 


(119) 


which  ARE  clearly  parameter  dependent!  Now  since  current  sensors  are  in  the  feedback 
path,  Kp„m  must  account  for  both  the  sensor  gain  and  the  gain  through  the  STPWM 

modulator  ( — ). 


D.  Extension  to  Two  Inverters  in  Parallel 

The  results  for  two  inverters  in  parallel  are  investigated  next.  Assuming  the  following 
parameters  for  the  system: 


L,  =  Lz  =  500gH 

(120) 

Ll  =  510 gH 

(121) 

Rl  =  4i2 

(122) 

and  that  the  single-unit  equivalent  connection  inductance  is  then 

L,  =  Y  =  250|iH 

(123) 

the  resulting  gains  are  (where  the  individual  inverter  control  gains  are  found  by  doubling 
the  result  since  we  have  two  in  parallel)  so  that 

Kpq,  =  K^qz  =  7.9373 

(124) 

K-q,  =  Kiqz  =  108,963 

(125) 

Kpdi  =  Kpdz  =  14.0506 

(126) 

Ki'd,  =  Kmz  =  86,863 

(127) 

To  derive  the  ACTUAL  gains  used  in  units  1  and  2,  the  PWM  and 

accounted  for 

sensor  gains  must  be 

T'  7.9373 

J^pql  -  J^pq2  -  ^ 

-^pwm  ^sensor 

(128) 

T.  108,963 

Kiql  -  Riq2  -  ^  ^ 

•^pwm  ■‘^sensor 

(129) 

,,  14.0506 

^pdl  ”■  ^pd2  “  ^  1^ 

•^pwm  -^sensor 

(130) 

r.  j.  86,863 

^idl  -  J^id2  -  ^  " 

^pwm  -^^sensor 

(131) 

With  (124)-(127)  substituted  into  (75),  the  following  actual  eigenvalues  are  achieved: 

^1,2  =  -7,943  +  j  12,445  (132) 
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^3,4  =  —5,258  +  j  6,642 

(133) 

^5,6  =  -7,238  ±  j  2,169 

(134) 

^7  =  -9,191 

(135) 

h  =  -18,899 

(136) 

Clearly,  there  is  a  strong  correspondence  between  the  desired  pole  locations  (106)-(107) 
and  those  realized.  Detailed  simulations  would  next  need  to  be  conducted  to  confirm  the 
dynamic  response. 


E.  Zero-Sequence  Control  Analysis 

The  zero-sequence  dynamics  for  two  inverters  supplying  the  R-L  load  are  given  by 


d  .  _  Vqj  -  Vo,2 

dt  ’  Lj  +  1^2 

(137) 

io,2  =  “io,i 

(138) 

Let  s  consider  that  only  Proportional  (P)  control  is  used  in  the  zero-sequence  loop,  so  that 

'^0,1  l^pwm  ^pO  (lo,l  lo,l  ) 

(139) 

Vq  2  —  Kpwm  l^pO  ( lo,2  “  *0,2  ) 

(140) 

With  io,i  =  io,2  =  0  and  (138),  it  follows  then  that 

Vo,i  -  Vo.2  =  -2Kp,^,mKpoio,i 

(141) 

which,  substituted  into  (137),  yields 

d  .  _  2  Kpwm  l^pO 

dt  ~  L,  +  L2 

(142) 

Equation  (142)  has  an  eigenvalue  at 

^  ^pwm  ^pO 

Lj  +  L2 

(143) 

and  thus  if  we  select  a  desired  eigenvalue  location 

the  required  proportional  gain  Kpo  is 

^  _  (L1+L2)  ^ 

(144) 

To  account  for  the  sensor  gain  for  the  zero-sequence  current  measurement,  (144)  would 
need  to  be  practically  restated  as 

^  _  (L]+L2) 

pO  0  XC  JC  ^des 

^pwm  ^sensor 

(145) 

The  A,<jes  should  be  at  least  on  decade  away  from  the  radian  switching  frequency  while 
sufficiently  far  in  the  left-half  plane  to  minimize  all  lower-frequency  circulating  currents. 
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F.  Extension  to  Five  Inverters  in  Parallel 

If  five  identical  units  are  placed  in  parallel,  the  expressions  for  the  qd-gains  become 


Kpq  =  5- 


(  ~  +  Ll  )  Kpq  ~  Rl 
l^pwm  t^sensor 


Kpd  =  5 


(y  +  LL)Kpd  -  Rl 
Rpwm  t^sensor 


K;„  =  5- 


(y  +Ll)K" 

Rpwm  Rsensor 


Kid  =  5 


(y +  Li,)Ki; 
Kpwm  Kgensor 


(146) 

(147) 

(148) 

(149) 


where  note,  the  connection  inductance  4  must  be  set  equal  to  the  coupling  inductance  of  a 
single  unit  divided  by  five.  The  resulting  gains  must  have  a  scale  factor  or  five  out  front 
for  reasons  documented  earlier  (in  particular  since  the  current  allocated  to  a  unit  is  one- 
fifth  that  of  an  equivalent  single  unit,  for  the  voltages  produced  to  be  the  same,  the  gains 
must  be  scaled  by  a  factor  of  five). 

To  address  the  design  of  the  zero-sequence  control  loops  for  five  inverters,  consider 
the  following  dynamic  equations 


“Vo,i  +  Lcyio.i  -  L(.  y  io,3  +  =  0 


(150) 

(151) 

(152) 

(153) 


where  an  identical  coupling  inductance  4  is  assumed  for  simplicity.  It  also  follows  that 
since  the  load  is  wye-connected,  one  zero-sequence  current  is  algebraically  related  to  the 
others.  If  we  choose  to  solve  for  io.j,  we  get 

■o..‘i  ~  ~io,i  ~  io.2  ~  'O-i  “  io,4  (164) 

If  we  assume  a  similar  proportional  control  for  each  zero-sequence  loop,  then  it  follows 
that 


vq.i  ~  Kpwm  Kpo  ( io,i  io,i ) 
Vo.2  “  KpwmKpQ(io^2  *0.2) 


(155) 

(156) 
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vo,3  =  Kp^Kpo(io,3  -  10,3) 

^^0,4  ~  ^pwm  ^pO  ( ^0,4  —  ^0,4  ) 
^0,5  ”  ^pwm  ^^pO  (  ^0,5  "  ^0,5  ) 


(157) 

(158) 

(159) 


Since  it  is  desired  that  Iq  ,  =  io2  =  io,3  =  io,4  =  io,5  =  0,  substituting  (154)-(159)  into  (150)-(153) 
and  placing  the  result  into  state-space  normal  form  yields 

0  0  0 


^0,1 

0  VmKpO  Q  Q 

*0,1 

*0,2 

Lc 

*0,2 

dt 

*0,3 

0  0  -Kpwn-^PO  0 

*0,3 

-^0,4 

Lc 

.*0,4„ 

l^wm  l^pO 


0  0 
The  characteristic  equation  for  (160)  is  simply 

KpwmKpO^^^  KpwmKpO  )  ( X.  +  )  =  0 

Lc  Lg  Lj.  Lf^ 

and  so  all  four  roots  are  co-located  at 

^1,2, 3.4  = 


l^pwm  l^pO 


(160) 


(161) 


(162) 


which  is  identical  to  (143)  with  L,  =L2.  Therefore,  the  zero-sequence  gain  formula  (144) 
would  continue  to  hold  for  five  inverters  in  parallel.  Repeating  that  equation  (assuming 
equal  combining/coupling  inductances) 


Kpo  = 


Lc  ^des 

Censor 


(163) 


where  the  impact  of  the  feedback  current  sensor  gain  has  been  incorporated. 
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IV.  INDUCTION  MACHINE  LOAD  ANALYSIS 

One  of  the  goals  of  this  research  is  to  facilitate  the  development  of  motor  controllers 
employing  multiple  paralleled  power  inverters.  The  following  analysis  will  extend  the  R-L 
load  results  to  an  induction  motor  drive  and,  in  a  later  section,  to  a  permanent-magnet 
synchronous  machine  drive.  The  stator  voltage  equations  of  the  three-phase  symmetrical 
induction  machine  represented  in  the  synchronous  reference  frame  are  given  by 


=  r 


'"qs 


*qs 


“b 


Vds 


_L1. 

(Ob  dt 


Vqs 


Vds 


(164) 

(165) 


where  Tj  is  the  per-phase  stator  resistance,  cOb  is  the  base  (rated)  electrical  angular 
frequency,  and  \|/  represents  flux  linkages  per  second.  The  stator  flux  linkages  per  second 
may  be  expressed  in  terms  of  the  stator  leakage  flux  and  the  air-gap  flux  by 


Vqs  =  (166) 

Vds  =  Xisils  +  (167) 

where  the  air-gap  flux  linkages  per  second  are  given  by 

Vmq  =  X,n(i®s  +  i'®)  (168) 

V.^d  =  X„(id*  +  i'^r)  (169) 


The  parameters  Xis  and  are  the  stator  leakage  and  magnetizing  reactances,  respectively 
(each  a  constant  evaluated  at  cob).  The  prime  in  (168)-(169)  denotes  that  a  rotor  variable  is 
a  referred  quantity  (that  is,  the  rotor  windings  have  been  mathematically  replaced  by 

equivalent  windings  with  the  same  number  of  turns  as  the  stator  windings).  The  referred 

rotor  flux  linkages  per  second  are  given  by 

Vq?  =  X,;i;«  +  V|/mq  (170) 

=  x,;i^  +  (171) 


where  X,',  is  the  referred  rotor  leakage  reactance  (a  constant).  Solving  (170)-(171)  for 
and  id^  substituting  the  results  into  (168)-(169)  and  collecting  the  and  \|/®d  terms  yields 


II 

X„x„ 

■qs  + 

Vqr 

x„  +  X,; 

X,; 

Vmd  = 

Xn,x,; 

ids  + 

Vdr 

Xn,  +  X,; 

X,; 

(172) 

(173) 


Finally,  substituting  (172)-(173)  into  (166)-(167)  gives 


Vqs  =  X"iq"3  +  \|/q 


(174) 


where  the  double-prime  quantities  are  given  by 

x„x,; 


X  =  Xu  + 


x„,  +  Xu 


Vq 


x„ 


x„  +  x„ 


% 


¥d  = 


x„ 


x„  +  X, 


- 


The  stator  voltage  equations  (164)-(165)  may  then  be  re-expressed  as 


1  d 


Vqs  =  +  —  (X  ids  +  ¥d)  + 


(176) 

(177) 

(178) 


(179) 


v|s  =  r,i|,  -  -^(Xi, 

®b 


+  ¥q  )  + 


_Li_ 

cOb  dt 


(X  +  ¥d) 


(180) 


The  current  control  dynamics  are  desired  to  be  very  fast.  As  a  result,  the  slowly-varying 
rotor  variables  may  be  assumed  to  be  constant  while  considering  the  stator  dynamics.  This 
implies  that 


d  'e  d 


=  — Vd?  =  0 


dt 


(181) 


which  observing  (177)-(178)  requires  that 


d  -  d  «  - 

=  ■JT'I'd  =  0 


(182) 


Furthermore,  assuming  the  drive  is  vector  controlled  ensures  that 
=  0 


¥d?  =  ¥* 


Substituting  (182)-(184)  into  (179)-(180)  yields 


to.. 


v^s  =  +  —  X  4  +  — 


tOe  X„ 


®b 


COb  Xm  +  Xir 


-¥*  + 


2Li_- 

(Ob  dt 


Aqs 


e  -e  “e,,".,  X  d  .. 
‘d-  * 


If  we  let 


L  =  — 


2^ 

tOb 


(183) 

(184) 

(185) 

(186) 

(187) 


and  merge  these  equations  with  the  qd-equations  for  a  single  inverter  with  coupling 
inductor  (81)-(82),  we  find 
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-  Vq  +  Ljj  iqs  +  Lx  COg  iqj  +  COg  L 


^e  ^m  'q*  -  "  d  ,  e  ^ 

'  Vdr  +  L  ^qs  “  ^ 


CDb  Xm  +  X|r 

— ’  V(|  +  Lx  ids  Lx  COg  iqs  +  Tg  i^jj  “  COg  L  iqs 

^  L-^iA  =  0 

If  we  next  observe  that  the  Pl-control  equations  are 

Vq  =  K^(iqt  -  iqs)  +  Kj,  X  | 

Vd  =  Kpd(i^  -  ids)  +  Kid^d 


dt 


—  X®  =  r 
dt 


*qs 


dt 


■”  ^ds  ^ds 


placing  (188)-(193)  in  state-space  normal  form  gives 


dt 


iqs 

-(Oe 

K" 

o' 

^qs 

ids 

COe 

-k;; 

0 

Ki^ 

ids 

Aq 

-1 

0 

0 

0 

Xq" 

XT  e 

0 

-1 

0 

0_ 

V  e 

^d 

Xd 

‘k;  o' 

1 

0  K;d 

1  -®eXn,V|f^* 

0 

1  0 

.‘ds* 

(Ob(Lx  +  L'')(X„,  +  X,;) 

0 

0  1  _ 

0 

where  in  this  case 

r,  +  K, 


Kpq  - 


Kpd  = 


pq 


Lx  +  L 

U  +  K,' 


pd 


Lx  +  L 


Kiq  = 


K;„ 


Lx  +  L 


Kid  = 


Kh 


L,  +  L 


(188) 


(189) 


(190) 

(191) 

(192) 

(193) 


(194) 


(195) 

(196) 

(197) 

(198) 


Note,  this  is  the  same  form  as  equations  (89)-(94)  where  the  stator  resistance  (rj  has 
replaced  the  load  resistance  and  the  stator  substransient  inductance  (l”)  has  replaced  the 
load  inductance.  Thus,  the  solution  may  progress  as  before:  use  the  Newton-Raphson  to 
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solve  for  the  double-primed  gains,  use  (195)-(198)  to  derive  the  primed  gains,  and  finally 
divide  those  quantities  by  Kp„m  Ksensor  to  establish  the  gains  for  the  single  aggregate  inverter. 
Those  gains  are  then  scaled  by  the  required  number  of  inverters  to  calculate  the  individual 
inverter  current  control  parameters. 


-  28  - 


IV.  PERMANENT-MAGNET  SYNCHRONOUS  MACHINE  LOAD  ANALYSIS 


The  derivation  of  the  control  equations  for  a  round-rotor  permanent-magnet 
synchronous  motor  (PMSM)  is  more  straight-forward.  The  PMSM  is  modeled  in  the  rotor 
reference  frame  by 

~  rsiqs  +  L55 iqj  +  (OfLsjijs  +  (199) 

ids  Tjs  ids  ~  l^ss  iqs  (200) 

These  equations  assume  negligible  cogging  torque  and  sinusoidal  back  emf.  The  parameter 
T,  is  the  per-phase  stator  resistance,  L^s  is  the  stator  self  inductance  (leakage  plus 
magnetizing),  and  ^  is  the  peak  flux  linkage  between  the  permanent  magnet  and  the  stator 
windings.  The  inverter  equations  in  the  rotor  reference  frame  are  similar  to  what  we  had 
previously  except  now  the  angle  in  (13)  is  the  rotor  electrical  angle,  superscript  ’e’  is 
replaced  by  superscript  ’r’,  and  cOe  is  replaced  by  the  rotor  electrical  angular  velocity  00^. 
Coupling  the  inverter  dynamics  in  the  rotor  reference  frame  (81)-(82)  to  the  machine 
equations  yields: 

~Vq  +  Lx  iqs  LxWfijjj  +  Tjiqj  (201) 

L55  iqs  +  03s  Ljs  ids  ^ 

-Vd  +  Lx -^4  -  Lx03si^s  +  Tsids  (202) 

Ljs  ids  ~  Lss  iqs  ~  9 

These  equations  may  then  be  combined  with  the  Pl-current  control  equations  in  the  rotor 
reference  frame 


Vq  ^pq  ^  iqs  iqs  1  i^iq  ^q 

"^d  kpd  (  ids  ids  )  kjd  Xd 

^  r  *  r*  • r 

^Xq  =  Iqs  -  Iqs 


dt 


Xd  =  ids  -  ids 


(203) 

(204) 

(205) 

(206) 


to  yield  the  state-space  description 
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^qs 

■-Kp;  -CO,  Ki;  o' 

iqs 

ids 

~  Kpd  0  IQd 

^ds 

-10  0  0 

X  ^ 

.  0  -1  0  0_ 

„  r 

1  Xd 

Kpq  0 

Lx  Ljs 

0  Kp'd 

0 

1  0 

0 

.  0  1  . 

0 

where  in  this  case 


Kpq  - 


^pq 
Lx  Ljs 


tc  +  K, 


"  Lx  +  L, 


•pd 


Ki„  = 


Ki. 


Lx  +  Ljs 


Ki,  = 


K. 


Lx  Lx; 


(207) 


(208) 

(209) 

(210) 
(211) 


Note,  once  again  this  is  the  same  form  as  equations  (89)-(94)  where  Rl  is  replaced  by 
and  Ll  is  replaced  by  Ljs-  Thus,  the  solution  may  progress  as  in  the  R-L  case:  use  the 
Newton-Raphson  routine  to  solve  for  the  double-prime  gains,  use  (208)-(211)  and  the 
equivalent  connection  inductance  to  solve  for  the  primed  gains,  and  divide  those  quantities 
by  KpvjTO  Ksensor  to  atrive  at  the  actual  gains  for  a  single  aggregate  inverter.  Those  gains  are 
then  scaled  by  the  required  number  of  paralleled  inverters. 
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VI.  VECTOR  CONTROL  IMPLEMENTATION 

Vector  control  is  the  state-of-the-art  means  of  controlling  AC  machines.  It  essentially 
gives  the  AC  machine  the  torque  characteristics  of  a  DC  commutator  machine  -  notably, 
near-instantaneous  torque  response  (as  quick  as  the  machine  currents  can  be  changed  and 
tracked).  The  following  equations  are  included  to  document  how  the  inverter  control 
development  presented  thus  far  can  be  extended  to  a  high-performance  servo  application. 
We  will  start  first  with  the  induction  machine. 

A.  Induction  Machine  Vector  Control 

The  induction  machine  vector  control  may  be  implemented  directly  (by  physically 
locating  the  position  of  the  rotor  flux  and  placing  the  synchronous  reference  frame  d-axis 
there)  or  indirectly  (by  using  the  induction  machine  governing  equations  to  establish  a 
necessary  and  sufficient  condition  to  guarantee  that  all  of  the  rotor  flux  is  directed  along 
the  synchronous  reference  frame  d-axis). 


1.  Direct  method  of  vector  control 

The  direct  method  of  vector  control  can  either  use  air-gap  flux  measurement  or  simply 
stator  voltage  and  current  measurements.  If  sense  coils  are  used  to  establish  the  air-gap 
flux  (by  integrating  the  resultant  coil  voltages),  then  we  can  establish  the  required 
synchronous  reference  frame  angle  by  the  following  sequence  of  calculations: 


^qs  — 

*as 

1  . 

2  . 

ids  = 

V3 

x;  , 

~  Xir  iqs 

¥qr  = 

^¥4q  • 

Am 

¥dr  = 

^Vmd  ■ 
Am 

-  x;,ii 

¥r,pk  = 

=  V(¥q?)^ 

+  (¥dr)^ 

(212) 

(213) 

(214) 

(215) 

(216) 

(217) 


Notably,  the  approach  has  issues  at  low-frequency  due  to  the  requirement  of  integrating  to 
find  the  air-gap  flux  linkage  components. 

An  alternative  scheme  does  not  require  special  air-gap  sensors,  but  only  simple 
current  and  voltage  measurements.  It  establishes  the  required  synchronous  reference  frame 
angle  by  first  calculating  the  stationary  reference  frame  voltages  and  currents  from  the 
physical  variables: 
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i  *  r=  i 
^qs  *as 

(218) 

ii 

1 

2  . 

V3 

(219) 

< 

II 

(220) 

11 

1 

2 

V3 

(221) 

The  stator  flux  linkages  per  second  (in  the  stationary  reference  frame)  are  then  found  by 
integrating  the  governing  equations: 

t 

¥qs  =  /(-rstObiqs  +  ®bV'3)d^  (222) 

0 

t 

¥ds  =  /(-rsCObils  +  ®bVd\)dC  (223) 

0 


This  again  imposes  a  DSP  difficulty  at  low  frequencies.  The  referred  rotor  flux  linkages 
are  then  determined  by 


¥qr  = 


5: 

x„ 


¥qs  ( Xjs 


¥dr 


x„ 


Vjs  -  (X,3- 


¥r.pk  =  V(¥qr)^  +  (¥dr)^ 

The  d-axis  is  placed  on  top  of 
assigning 


x^ 


(224) 

(225) 

(226) 


the  rotor  d-flux  in  the  synchronous  reference  frame  by 


=  tan  ' 


(227) 


2.  Indirect  method  of  vector  control 

Alternatively,  the  desired  synchronous  reference  frame  angle  may  be  derived  by  only 
measuring  the  rotor  electrical  angular  position,  where  this  relates  to  the  mechanical  angle 
by 


=  yerm  (228) 

where  P  is  the  number  of  poles.  Then  Ge  is  established  from  a  necessary  and  sufficient 
condition 


(Obr^X, 


•  e. 

m  *qs 


Xrr¥** 


dC 


(229) 


where  vt/jf  is  the  desired  d-axis  rotor  flux  and  i®*  is  the  desired  q-axis  stator  current,  both  in 
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the  synchronous  reference  frame.  If  the  desired  rotor  flux  magnitude  is  kept  constant,  then 
it  follows  that 


Vdf  =  x^ii; 


(230) 


and  (229)  simplifies  to 


COfeTri, 


qs 


Xrril. 


dC 


(231) 


Thus  it  requires  only  position  feedback  information  and  feedforward  of  the  commanded 
stator  currents  in  the  synchronous  reference  frame. 


B.  Permanent-Magnet  Synchronous  Machine 

Field-oriented  vector  control  for  a  PMSM  is  achieved  through  self-synchronization 
(feeding  back  6^  and  setting  it  equal  to  e*)  and  setting  equal  to  zero.  The  commanded 
q-axis  stator  current  will  then  be  directly  proportional  to  the  developed  electromagnetic 
torque. 


-  33  - 


VII.  OUTER  SPEED  LOOP  DESIGN 


A  perspective  outer  speed  control  loop  for  a  vector-controlled  induction  machine  or 
PMSM  may  take  the  form  illustrated  in  Figure  7.  That  is  under  vector  control,  the 
developed  electromagnetic  torque  is  related  to  the  commanded  q-stator  current  for  an 
induction  machine  by 


Te  = 


3PXj^.e* 

4(0bx; 


(232) 


and  for  a  round-rotor  PMSM  by 


T  = 


4  ‘qs 


(233) 


Figure  7.  Closed-loop  speed  control  system 


By  inspection  of  (232)-(233),  the  appropriate  value  for  Kt  in  Figure  7  may  be  found.  The 
closed-loop  transfer  function  for  the  control  loop  of  Figure  7  is  given  by 


Gcl(s) 


_ K,Kp,s  +  K,Ki, _ 

Jtot  +  (  Btot  +  Xt  Kps )  s  +  K,  Kjj 


(234) 


where  Jjot  represents  the  combined  inertia  of  motor  and  load,  Bjot  is  the  combined  viscous 
friction  or  motor  and  load,  and  Kps  and  Kis  are  the  Pl-control  gains.  By  specifying  the 
dynamics  of  the  speed  loop  (through  a  pair  of  desired  pole  locations),  the  two  gains  may 
be  uniquely  determined.  The  pole  locations  are  best  selected  at  least  one  decade  away 
from  the  current  control  pole  locations  to  avoid  unwanted  interaction.  The  zero  in  (234) 
may  impact  the  response  unfavorably  by  introducing  more  overshoot  than  desired.  One 
potential  solution  is  to  insert  a  pre-filter  after  the  commanded  speed  signal.  For  a  pre-filter 
choice  of 
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Gpf(s) 


1 


K, 


•ps 


K. 


+  1 


(235) 


the  zero  is  canceled.  Obviously,  detailed  simulation  should  be  used  in  conjunction  with 
these  design  suggestions  to  investigate  interaction  and  non-idealities  in  full. 
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VIII.  CONCLUSION 

The  research  effort  documented  in  this  technical  report  provides  analysis/synthesis 
tools  for  choosing  feedback  gains  for  a  set  of  paralleled  inverters.  The  local  control  is  set 
forth  in  a  synchronous  reference  frame  and  also  addresses  the  question  of  circulating  zero- 
sequence  currents.  The  procedure  for  a  generic  R-L  load  was  extended  to  both  an 
induction  motor  load  and  a  permanent-magnet  synchronous  machine  load.  Interestingly,  a 
similar  set  of  equations  was  found  to  be  applicable  in  each  case.  The  resulting  inner 
current  control  design  was  then  mapped  outward,  in  the  case  of  the  machine  loads,  to 
facilitate  the  design  of  a  speed-control  drive.  To  make  the  analysis  comprehensive, 
equations  for  the  vector  control  implementation  options  were  presented  and  the  resultant 
simple  Pl-control  loop  analyzed. 

In  order  to  validate  the  control  procedures  outlined  here,  the  dynamic  equations  for 
the  R-L  load,  induction  motor,  and  PMSM,  together  with  the  inverters,  must  be  simulated 
in  detail  using  a  program  such  as  ACSL.  Such  simulations  could  be  structured  to  model 
slight  differences  between  the  units  and  thus  study  the  efficacy  of  the  circulating  current 
control.  Also,  the  interaction  between  sensor/filter  dynamics  and  the  current  control 
dynamics  may  be  more  accurately  studied  and  a  more  informed  decision  made  about  the 
suitability  of  the  current  control  bandwidth.  Finally,  simulation  would  allow  the  study  of 
the  vector  control  implementations  and  any  sort  of  speed  control  loop.  The  preceding 
research  effort  must  next  be  applied  to  a  hardware  testbed,  where  the  individual 
programming  the  DSP  control  must  work  closely  with  the  modeler/simulator  to  accurately 
test  each  stage  of  development. 

Ultimately,  the  above  approach  could  be  extended  to  a  propulsion-type  application  of 
a  high-bandwidth  low-power  inverter  module  (consisting  of  paralleled  units)  operating  in 
parallel  with  a  low-bandwidth  (slower  switching)  high-power  inverter  module.  The 
efficacy  of  using  the  high-bandwidth  unit  to  improve  the  fidelity  of  the  output  waveforms 
could  then  be  evaluated. 
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